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Abstract 

We study the hermitian one matrix model with semi-classical potential. This 
is a general unitary invariant random matrix ensemble in which the potential has 
a derivative that is a rational function and the measure is supported on some col¬ 
lection of disjoint closed intervals. Such models have attracted much interest both 
due to their physical applications and relations to integrable systems. An object of 
central interest in random matrix theory is the correlation kernel, as this encodes 
the eigenvalue correlation functions. In recent years many results have been ob¬ 
tained proving that the correlation kernel near special points in the spectrum can 
be expressed in terms of Painleve transcendents and their associated Riemann- 
Hilbert problems. In the present work we build on this success by proposing a 
model problem that is general enough to describe the limiting kernel at any point 
in the spectrum. In the most general situation this would include cases of loga¬ 
rithmic singularities and essential singularities in the weight colliding with soft or 
hard edges, the bulk of the spectrum or even births of a cut. 


1 


Contents 


1 Introduction 

1.1 Motivations. 

1.2 Statement of results . 

1.2.1 The canonical RH problem for . . . 

1.2.2 Definitions for the double scaling limit 

1.2.3 Main theorems. 



2 Properties of the canonical model problem 

2.1 Opening the lens. 

2.2 Examples. 

2.2.1 Relation to the Airy model problem .... 

2.2.2 Relation to the Bessel model problem ... 

2.2.3 Relation to Bessel in the bulk. 

2.2.4 Relation to the general Painleve II equation 

2.2.5 Relation to GUE at a birth of a cut .... 

2.3 Solvability of 4> RH problem: Proof of theorem 11.241 


Q 

13 

15 

16 
16 

U 

18 

18 

18 


3 


Asymptotic analysis of the RH problem for orthogonal polynomials 

3.1 The RH problem for orthogonal polynomials. 

3.2 Transformation to constant jumps . 

3.3 Opening the lens. 

3.3.1 Definitions of contours. 

3.4 Normalisation at infinity. 

3.5 Final transformation. 

3.6 Proof that R is a small norm RH problem. 


23 

23 

24 

25 
25 
27 
27 
3C 


4 Construction of global and local parametrices 

4.1 Construction of . 

4.2 Asymptotic behaviour of ^ and near singular points 

4.3 Construction of . 


1 ^ 

32 

38 


5 Asymptotics for Kn{x,y)'- Proof of theorem 11.261 


41 


6 Example applications 

6.1 The Painleve H kernel in the bulk 

6.2 A new kernel at a hard edge . . . 


43 

44 


2 














































1 Introduction 


Let X be a finite collection of closed pairwise disjoint intervals in M and 1/ : M —>• M 
a function whose derivative is rational. Then define a semi-classical hermitian matrix 
model as a random n x n hermitian matrix given by M = C/ldiag(xi,... ,Xn)U where 
[/ is a Haar distributed unitary matrix and the eigenvalues Xj are distributed according 
to the probability measure, 

1 ” 

(1.1) 

i=l 

Here xi is the indicator function on I, A{x) is the Vandermonde determinant, the 
function V is known as the potential and we often refer to w{x) := as the 

weight. The normalisation constant 

~ n 

Z„= / A{xf\[e-^^^^^hx{xi)dxi. (1.2) 

i=i 

is known as the partition function. Note that implicit in this definition is that V is 
such that Zn exists. 

The eigenvalues distributed according to o form a determinantal point process. 
The correlation functions for this process can therefore all be expressed in terms of a 
correlation kernel. In turn the correlation kernel is expressed using orthogonal polyno¬ 
mials with respect to the weight w{x) on I. More precisely, let pj, j = 0,1,... be the 
family of monic polynomials of degree j characterised by the relations 

j Pj{x)pm{x)w{x)dx = hjdjm- (1-3) 

The correlation kernel can now be written as 

Kn{x,y) = (p„(x)p„_i(y) - pn{y)pn-i{x)) . (1.4) 

X y 

The asymptotic analysis of the correlation kernel therefore reduces to the asymptotic 
analysis of the orthogonal polynomials pn- Indeed, the use of the term semi-classical 
to describe such an ensemble arises from the fact that their associated orthogonal 
polynomials are of semi-classical type. 

In the present work we will obtain n —)• oo asymptotics for the eigenvalue correlation 
kernel of such a model under very general conditions. In particular we obtain the 
asymptotic behaviour of the kernel at points in which edges of X and singularities in 
the potential approach the support of the eigenvalues as n —)• oo in a way that ensures 
a non-trivial double scaling limit. We will express the asymptotics for the kernel in 
all such cases in terms of a general model Riemann-Hilbert problem which generalises 
nearly all other model problems previously introduced in the context of random matrix 
theory. Some examples of the scenarios not addressed in the analysis here include cases 
when the potential has isolated discontinuities, which have been studied in [T] and 
lead to confluent hyper geometric kernels and also the interpolating kernel found in the 
context of the birth-of-a-cut [2]. Finally, for reasons of space we restrict ourselves to 
the case of logarithmic singularities in V with positive coefficients. 
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1.1 Motivations 


We begin this section with a historical review of the physical applications of random 
matrix theory, after which we turn to the mathematical motivations for the current 
study. 

The one hermitian matrix model in which V is polynomial and I = R was the first 
random matrix model to be extensively studied. These studies culminated in signihcant 
applications to problems in particle physics [3]. Further applications quickly followed, 
with many requiring some modification of the model to allow V to be rational or to 
allow / 7^ M. To be more precise it was found that for applications of random matrix 
theory to QCD [3], quantum transport problems [HIEI and string theory [3, the models 
necessarily included logarithmic divergences in V. 

Since that time there has been an explosion in the number of applications of random 
matrix theory and many of the proposed models fall into the class of semi-classical 
hermitian matrix models. This is especially true in the context of quantum transport 
in which recent models incorporating poles in to the potential have been shown to have 
physical relevance. There the observable of interest is the Wigner-Smith time-delay 
matrix Q, whose eigenvalues tj are related to the “inverse delay times”, The 

joint probability density for the was hrst obtained in [5l [6] where it was shown it 
took the form. 
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i=i 


(1.5) 


where j3 depends on the symmetries of the system, with /3 = 2 a common case. Since 
many observables may be expressed in terms of Q, the problem of computing expecta¬ 
tion values with respect to the above measure is relevant. Let us highlight the recent 
work [8] and [9] in which physically important observables were expressed in terms of 
Tr Q and Tr Q^. The associated moment generating function for these observables is the 
partition function for a semi-classical hermitian matrix model in which the potential 
contains a singularity of order k = 1 and k = 2 respectively. 

Random matrix theory has also found applications in integrable quantum field the¬ 
ory at finite temperature. Here again it was found that models of semi-classical type 
were relevant [lOj . 

Finally, another important application of random matrix theory is in the held of 
analytic number theory. Here random matrix theory serves as a source of models for 
the zeros of various C-hke functions. An example of this relation was given recently by 
[?] in which the statistics of large deviations of \C,{l/2 + ix)\ in an interval on the critical 
line was found to be related to the statistics of large deviations of | det([/ —e*®I)|, where 
U is a n X n Harr-distributed unitary matrix. 

Of central importance in the work m was a conjecture for the asymptotics of 
certain expectation values using a measure of the form (11.11) but where the measure is 
supported on a circle in the complex plane rather than the real line. Such ensembles are 
known as Circular Unitary Ensembles (CUE). The CUE is very closely related to the 
semi-classical hermitian matrix model and we expect the same set of model problems 
as introduced here to appear in the CUE context. To address the conjecture made in 
m would require analysing model problems with an arbitrary number of logarithmic 
singularities, which is a property of the model problems introduced in this paper. This 
is something we hope to pursue in the near future. 
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Another work related to the C function is [T^] . There the motivating interest was in 
a random matrix model describing a certain observable with relevance to the Riemann 
hypothesis. The authors of |12j studied a semi-classical hermitian model which contains 
a simple and second order pole . 

We now turn our attention to motivations arising from Riemann-Hilbert analysis. 
The use of Riemann-Hilbert problems and the associated machinery of the Deift-Zhou 
steepest descent has enjoyed considerable success in random matrix theory over the 
last fifteen years. R has been able to give rigorous proofs of many conjectures made in 
the physics literature. We highlight especially the strong universality results that have 
been obtained using these methods. For example it has established universality of the 
Sine, Airy, Bessel and numerous Painleve kernels, for a broad class of potentials. 

The industry of applying the Deift-Zhou steepest descent to the RH problems arising 
in the context of random matrix theory is still strong with numerous results in the last 
few years da [mini [IS]- These results in particular have begun to address the types 
of new critical behaviour that one finds in the semi-classical hermitian matrix model. 
The current work was motivated by the desire to continue this program. In doing so 
we obtain proofs for many small conjectures in the Riemann-Hilbert analysis literature. 
For instance our main results prove the conjectures made in [2| and m concerning the 
most general behaviour at the birth of a cut and at a critical edge point. 

Finally let us note that there exists an extensive literature dzj on an alternative 
method for computing asymptotics of random matrix integrals, known as topological 
recursion. Topological recursion gives a uniform structure to the asymptotic expansion 
of observables in random matrix theory. Furthermore, the recursion needs only the 
data of the spectral curve together with a geometric object known as the Bergmann- 
kernel. This work was partly motivated by a desire to better understand the role of 
these objects in the Riemann-Hilbert approach. In particular our results show that 
the spectral curve directly determines the double scaling limit in all cases, without the 
need for the kind of smart approximations made in [laiaia. The Bergmann-kernel is 
a bit more mysterious in the current formulation but this information likely resides in 
the global parametrix together with the RH problem for the orthogonal polynomials. 

1.2 Statement of results 

Our results are the following: 

1. We introduce a general model RH problem which we refer to as the canonical 
model problem. This model problem generalises many of the model problems 
already used in the random matrix theory literature, such as the Airy, Bessel, 
Painleve I and H model problems. The canonical model problem has jump ma¬ 
trices independent of the parameters of the system and therefore possesses an 
associated system of linear ODEs which corresponds, when it exists, to the Lax 
pair of a Painleve type equation. 

2. We prove the existence of a solution to the canonical model problem for a certain 
range of model parameters. 

3. We prove, given a suitable definition of the double scaling limit, that the double 
scaled correlation kernel of the eigenvalues may always be expressed in terms of 
a solution to the canonical model problem. 
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Figure 1: An example of a system of contours in a canonical model problem. The thick 
lines are the intervals in I. The thin dashed lines are the intervals in IJfegB ^b\I- The 
black dots, one of which is labelled b, are locations of points in B and the contours 
ending at these points are the F^ contours. The shaded grey regions denoted by 
correspond to the potential location of sectors in which the asymptotic behaviour of $ 
differs. 

We will state our results more precisely in Theorem 11.241 and Theorem 11.261 This 
will require that we first introduce some more notation and definitions in order to define 
the canonical model problem and the double scaling limit. 

1.2.1 The canonical RH problem for 

Let us begin by motivating the following definitions. The model problem will arise as 
the local behaviour of a RH problem for orthogonal polynomials associated with the 
weight in (II.ip . Naturally therefore, the model problem must contain the features from 
the weight in (11.11) . In particular it must contain a collection of closed disjoint intervals 
I which encodes the local behaviour of I together with a set B which represents the 
singular points of V and a set of vectors {rfe} giving the coefficients of the singular 
terms in V. With this in mind we now give some definitions. 

Definition 1.1. Let R be a finite set of points 6 G C together with their complex 
conjugates. To each b £ B we associate a vector fj, of length db- 

Because the singular points in b may contain logarithmic branch points we must 
also define the jump contours ending at each b. 

Definition 1.2. Given a point 6 G C, define the contour, 

F{, := (—oo. Re b) U (Re b + f[0, Im b]) (1.6) 

with the orientation of the contour directed away from —oo. 
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Definition 1.3. Let [z — 6]^“ denote the function (z — 6)^" in which the branch cut is 
taken on Lf, and [z — 6]^" is positive for z — b G M'*'. 

We need a function on to describe the jumps due to the logarithmic di¬ 

vergences. We give this in the following definition together with an extension of the 
definition into the complex plane. For convenience define dfe := Tb^ G [0,oo) for each 
b G B. Then introduce, 

Definition 1.4. 


ar{z) 


Ylb&BXrd^)ab, zG\Jf,^BTb 

'Eb&BXri,(^e{z))ab, z G C\[j^^j^{Reb + iR). 


Finally we define the very useful function, 


Definition 1.5. 6{z) := ±1 for ±lmz > 0. 

We are now in a position to define the relevant class of model problems. 


(1.7) 


Definition 1.6. Let / be a collection of pairwise disjoint closed, possibly infinite, 
intervals. Let be a vector whose length we will specify later. We define a canonical 
model problem for a function B, {Tb},Too) to be the following RH problem. Note 

that in the following we suppress any unnecessary arguments of 4>. 


RH problem for 

(a) $ : C \ / \ UftgsFfe —>■ is analytic in 

(b) The jump matrix j(z) := 

Xi{z) \ 

J\Z) 1 g g-27riOr(2) j ’ 


z. See Figured) 

has the following form, 

zGlu[jTb. 

b&B 


( 1 . 8 ) 


(c) To describe the asymptotic behaviour of as z ^ oo we need the function. 


k{z;n) 




( 1 

y-6i(z)e2™r0(2) 



where is a sector in C and 9 is defined in 
Define a polynomial of degree k by, 
k 


for z G ^ 
for z G C \ D. 
Definition 11.51 


Pk{,z) ■— ^ ^ BqojZ^ . 
1=0 


(1.9) 


( 1 . 10 ) 


We then have three possibilities as z —>■ oo, 

(i) we say is exterior-type of order-A: with behaviour, 

^^{z) = (/ + 0(z-^)) 
where c G N, fc G N and dtot = Ylb&B 
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(ii) we say is edge-type of order-A; with behaviour, 

where N = -^{I -\- iai), k G {—1} U 2N° and Too,fc+i G for k G 2N° and 

Too,k+i £ for k = —1. We also introduced the sector = { 2 : : TT > 

I arg^l > TT —7r/(2/i; + 3)} if A: G 2N° and O^®'^®®) = { 2 ; : vr > | argzj > tt —tt/S} 
when A; = — 1. 

(hi) we say is interior-type of order-A; with behaviour, 


(1.13) 

where k G N*^, roo, 2 fc+i £ = {z : \ argz\ < 7r/(4A: + 2)}, = 

{z : TT > \ arg^l > vr — 7r/(4A: + 2)}, 9 is dehned in Dehnition 11.51 and 


Q{z) 



Imz > 0 
Imz < 0. 


(1.14) 


In the above the principal branches of 2 ;^/^, 2 and log 2 are taken analytic off 
(— 00 , 0 ] and positive for 2 > 0. 

(d) As 2 —)• 6 G B, 


<h( 2 ) = 0(1} exp 
As z ^ a € dl \ B 


di-l 


i=i 


[z-b] 


afecra 


(1.15) 


^ , , fe27rii°g(^ “)o'+ if a is a right edge of an interval, , 

$( 2 ) = $a,o(^) X < _^log(a-z)a+ -r ■ if. i f i 

I e 27ri j if a IS a left edge of an interval. 


Here <I)a^o is analytic with respect to 2 at a. 


Remark 1.7. The above model problem has jump matrices which are constant with 
respect to all parameters appearing in the problem. This has the consequence that Lax 
matrices may be constructed from the logarithmic derivatives, cl4><l>“^, with respect to 
any parameter appearing in the problem. The compatibility of these derivatives then 
leads to a system of ODEs which will possess the Painleve property. Hence the model 
problem is a Riemann-Hilbert problem for ODEs of Painleve type. 


1.2.2 Definitions for the double scaling limit 

Definition 1.8. Recall the definition of the model (11.11) . Let A := dX. This set 
contains the endpoints of the intervals comprising X. 





Definition 1.9. Given a singular point z of ^ we write the principal part of the 
Laurent expansion of ^ at whose inner radius is zero, as 

d^—l 

tz,j{z - (1.17) 

j=0 


Let B be the set of all singular points of Since V'{z) is a rational function we may 
parameterise it as, 


V(z) : — V[eg{z) + l^ing(2!) + Lbi.(2;), 


where. 


fioo + 1 


yregiz) 


■Z^ 

00,J'^ ? 




hsing(2) Y^ 


beB 


db-1 




i=i 


E 2 

-ab log \z-b\. 
n 


bets 


(1.18) 


(1.19) 

( 1 . 20 ) 
( 1 . 21 ) 


Note that the fact that the probability measure for the model must be real implies b G 
B b* G B. For convenience we define some shorthand notation; w{z) := 

Wregiz) := rcsing(2) := Whr{z) := We restrict ourselves 

to the case that ab is a positive constant. Finally we define l^eg such that. 


lim 

x^+oo log(x2 + 1) 


= + 00 , 


( 1 . 22 ) 


and that l^ing is such that the integral exists. 

Definition 1.10. The measure dfi{x) is defined as the equilibrium measure which 
minimizes 

H^^) = jj \ y^ d^x{x)dp,{y) + j Vreg{y)dn{y), (1.23) 

among all Borel probability measures /r on I. The equilibrium measure can be written 
in terms of a density p; dp{x) = p{x)dx. Define S := supp p. 


Remark 1.11. The equilibrium measure satisfies the following inequalities. 


2 J log |x - y\p{y)dy -V{x) =£, x €S, 

2 J log \x — y\p{y)dy — V{x) < £, x Gl\S. 


(1.24) 

(1.25) 


Definition 1.12. The inequality (|1.25l) is taken to be strict for x £ I\S outside a 
finite number of isolated points. We define £ C X \ 5 to be the set of such isolated 
exterior points and 77 := 5 U T. Finally let p := sup 77. 
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Remark 1.13. Since the equilibrium measure is defined using l^eg the behaviour of p 
is well known (201 ED and can be expressed as, 

P{x) = ^y(x)+, (1.26) 

where we have introduced the spectral curve, 

y{z) := h{z)^/R{z). (1.27) 

Here h is a real polynomial and i? is a rational function with only simple poles. The 
set of zeros and poles of R coincides with dS and the function \jR{z) is defined with 
branch cuts coinciding with S. In (I1.26P we have taken the value of sjR{z) on the 
positive side of the cut. The behaviour of p is best discussed in terms of three distinct 
types of points: 

• An edge point x of order k € {—1} U 2N*^ is a point such that x £ dS and 

y{z) = 0{{z — x)^^^) asz^x. (1-28) 

A singular edge is one whose order k > 0, a hard-edge is one for which k = —1 
and a soft-edge is one with k = 0. 

• An interior point x of order /c > 0 is a point such that x £ S and 

y{z) = 0{{z — x)"^^) asz^x. (1-29) 

A singular interior point is one whose order A: > 0. 

• An exterior point x of order A: > 1 is a point such that x £ £, i.e. a point at 
which the ()1.24|1 holds outside of S, and 

y{z) = 0{{z — x)‘^^~^) as z—>■ x. (1.30) 

Such points are also known as “a birth of a cut”. 

Definition 1.14. Let 7Z be the set of zeros and poles of R and R be the set of zeros 
of h. 

Definition 1.15. Define the “^r-function” 

g{z) := j log( 2 : - x)dp{x), (1.31) 

where the principal branch of the logarithm is taken, meaning g is analytic on C \ 
(—oo,p]. The ^'-function has a number of properties we will make use of, in particular 


9 +{x) + g-{x) -v{x) -£ = 0 

X £ , 

(1.32) 

9 +{x) + g-{x) -V{x) -£<Q 

X £l\J, 

(1.33) 

rp 

9+ix) - g-ix) = 2Tri / p{x)dx. 

J X 


(1.34) 

Let us also define, 



^( 2 ;) = -^/ h{s)^/R{s)ds 

for z G C \ (—oo,p]. 

(1.35) 
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where the integration path does not cross the real axis. By (| 1.341) we have, 


9+{x) -g-{x) = 2C+(x), 

(1.36) 

and usine (11.321). we obtain the identitv 


= 2g{z) - Viz) - L 

(1.37) 

The jumps of ^ follow from those of 5 , 


^+(x) +C-ix) = 0 X G J, 

(1.38) 

^+{x) + C-ix) < 0 xGl\J, 

(1.39) 

rp 

^+(x) — C-ix) = 27ri / p[x)dx. 

J X 

(1.40) 

It is imoortant to note that (11.401) has the conseouence that. 


^+(x) — ^_(x) = 27rie(x), x € M \ 

(1.41) 


where e : M —)• [ 0 , 1 ] is constant on each connected component of M \ and is known 
as the “filling fraction”. 

Assumptions 1.16. The sets A, B etc. are n dependent. For an any hnite set X of 
n dependent points we assume that the limit. 


A* := lim A (1.42) 

n—)-oo 

exists. For some x G X we also use the notation x* := lim„^oo x. We extend these 
definitions to intervals and collections of intervals by defining [a, b]^ := [a*, 6 *]. We also 
assume that A* C ^7* i.e all singular points that are due to the weight collide with 
the support of the eigenvalues as n —oo. 

The singular behaviour of the weight in (II.ip is not the only way in which singular 
points of the spectrum arise. Indeed the most common way is when a zero of h ap¬ 
proaches J. We therefore need to keep track of these zeros of h which collide with J 
as n — 7 > oo. We do this via the following definition. 

Definition 1.17. Let us define % to be the largest subset of % such that "H* C J^. 
It will be useful to collect all the points that cause singular behaviour, besides zeros of 
7^, together in a single set, 

V ■.= AyjByj'H. (1.43) 

Definition 1.18. For x £ P U A let x* := lim^^oo x. Given Assumption 1 1.16) we have 
X* € and therefore that x* can be classified according to the scheme in Remark 1 1.1 31 
We say x scales appropriately as n —)> oo if |x — x*| = where is defined 

by, 

if X* is a edge point of order k, 

if X* is an interior point of order k, (1-44) 

if X* is an exterior point of order k. 

Furthermore, given b G B we say the parameters 4 scale appropriately if lim^^oo 4 
exists. 




2k+3^ 

1 

2 fc+l’ 

J_ 

2fc’ 
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Assumptions 1.19. We assume all points in P U 7^ scale appropriately. 

Definition 1.20. Given a set A C VUTZ and x* £ P* U7^* we define X\x, := {z £ A : 
2 * = X*}. For a collection of intervals I we define I\x^ to be the collection of intervals 
such that i9(/|a;*) = {dl)\xt.- As it stands, the definition of I\x, is ambiguous as it only 
specifies the endpoints of the intervals. We take care of this by further requiring that if 
X* is an exterior point, then I\x^ does not have intervals extending to ±oo, if x* is an 
interior point then I\x^ has intervals extending to both ±oo and if x* is an edge point 
then I\x, extends to only — oo. 

Definition 1.21. For x* £ P* U 7^*, define f{z) := n^^*{z — x*). The definition of / 
extends to intervals by defining /([xo,xi]) := [/(xq), /(xi)] with similar definitions for 
intervals with open ends. 


Definition 1.22. We define the scaling limit of y at x* by. 


y(C) := lim ^''*y(x*+n 

n^oo 


Define also, 

1(C) := 


K 


yiz)dz, 



X* is a exterior or edge point, 
X* is an interior point. 


(1.45) 


(1.46) 


where p := lim^^oo f{p') and p' := supT^-lj,^ if supT^-lj,^ exists else p' = x*. 

Definition 1.23. We define the model problem atx^, to be the model problem 

of the same type and order as the point x*. The vector Too is determined using the 
scaling limit of y at x*. We choose Tqo in the following way, 


(i) If X* is an interior point of order k we have for large Q, 

2k+l 

7=0 


(1.47) 


where Ej are constants. We set Too,o = Eq + (2k + l)(nf^+(x*) + 7ra+^) and for 
j > 0, Tooj = Ej. 


(ii) If X* is an edge point of order k we have for large C; 

k-\-l 

CHiC) = + ( 1 - 48 ) 

7=0 

where Ej are constants. We set Tooj = Ej. 


(iii) If X* is an exterior point of order k we have for large C, 


1 


2k 


m = -clog c - ^ E ^iC^'+ o(c^), 

7=0 

where Ej are constants. We set TooJ = Ej. 

The remaining data for the model problem is, 

1 = lim f(I\x,), 
n^oo 

B= lim f(B\x,), 

n^oo 

n = lim n^^*tf-i(h). 


(1.49) 
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(1.50) 

(1.51) 

(1.52) 




1.2.3 Main theorems 

Theorem 1.24. A solution of the canonical model problem for exists in the edge 
and interior ease for any real vectors Tb subject to the constraint that i/ 6 £ n M there 
exists a constant <5 £ M such that, 


( 

\z-b 


db 


> 5 , 


Vz £ I. 


(1.53) 


When $ is of exterior type of order k a solution exists if the above constraint is 
satisfied and 


J e (1-54) 

exists. 

Remark 1.25. In (jl.53l) recall that db is the order of the pole in V appearing at b. 
The condition (I1.53P ensures that if an essential singularity of the weight is approached 
along a contour in I, <I> remains bounded. The extra condition on <I> in the exterior case 
is related to the fact that in this case may be constructed explicitly from orthogonal 
polynomials. 


The second theorem concerns the behaviour of the kernel near a point x* £ P* U7^*. 

Theorem 1.26. Consider a semi-classical matrix model with ab > 0 for all b ^ B. 
Given a point x* £ P* U 7?.* define the analytic functions fii : 'E. \ I \ Ub^B^b C by, 


\Mz)) 



(1.55) 


where we have used the model problem at x* defined in Definition \1.23[ We also define 
the ^-kernel. 


a') ■= r~'^har(u)+ar(v)) 4>l{v)4‘2{u) 

^ ^ ■ 27ri(u - v) 


(1.56) 


In the double scaling limit where n —oo such that all parameters scale appropriately 
we have that. 


lim n ^^*Kn(x,t + n x*+n = K‘*’(m, u), 

n—)-co ^ ^ 


(1.57) 


foru,v £ lim„^oo/(2iUJ- 


2 Properties of the canonical model problem $ 

2.1 Opening the lens 

In order to compare to known model problems it is useful to transform the model 
problem by opening the lens. To this end we introduce some lens-like contours. 

Definition 2.1. If <I> is of interior type let L := {—oo,pi) U {p 2 ,oo) C I for pi < 
inf ReR < supReR < p 2 . If d) is of edge type let L := (—oo,pi) C I withpi < inf ReR. 
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Exterior Type 



Figure 2: Examples of jump contours for the canonical model problem with open 
lenses. The thick lines are the intervals in I. The thin dashed lines are the intervals in 
I U UbgB Tb \ L. The thick dashed line are the intervals in L. We have also labelled the 
lens contours and the lens regions The black dots are locations of points in B 
and the contours ending at these points are the T;, contours. 


The set L will dehne the parts of the real line enclosed in the lens. Note that in the 
above definition we required that the lens not include any points from B. 

Definition 2.2. If $ is of interior type let s||_ be rays inside from pi to 

—oo such that ±Im 2 > 0 on each ray respectively and orientated away from infinity. 
Let be rays inside from p 2 to oo such that ±Imz > 0 on each ray 

respectively and orientated towards infinity. Finally let := s||_ u 

Definition 2.3. If <I> is of edge type let be rays inside from pi to —oo such 

that ±Imz > 0 on each ray respectively and with orientation away from infinity. 

Definition 2.4. We define to be the region of C bounded by LU such that the 
+ side and — side of L is contained in and respectively. 

We now make the transformation, 

^k{z\I,B,{n},Too) := ^kiz\I,B,{Tb},Too)K{z,n), ( 2 . 1 ) 

where Q = U This leads to the RH problem, 

RH problem for 

(a) $ : C \ / \ \ \ —>■ is analytic in z (see Figure [2|). 
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( 2 . 2 ) 


(b) The jump matrix j{z) := ‘h_( 2 :) ^<h_|_( 2 :) has the following form, 


Jiz) = 

Ji^) = 

j{z) = 


0 1 
-1 0 
g2*7rar(^) 

0 g-2i7ror(2:) I ’ 


1 0 

g±2i7rar(2) ]_ 


z G T, 

^g/u (^r;,\L, 


b&B 


z G 


(2.3) 

(2.4) 


(c) As z ^ 00 , 

(i) If 4) is of exterior-type, 

$(Z) = [I + 0{Z-^)) 2(«t°t+^)<^3g-^T2fc(^)^3, 

(ii) If <h is of edge-type, 

$(^) = (7 + 0 ( 2 "^)) 

(hi) If 4* is of interior-type, 

$(z) = (7 + 0(2"^)) e-5rTT^2fe+i(^y3g(_2)_ 

Note that we have made use of the fact that = e°‘'^^Q{z) 

(d) As z —)• 6 G B, 

d(, —1 


(2.5) 


( 2 . 6 ) 


(2.7) 


<h(z) = 0(1) exp 
As z ^ a & dl \ B 


i=i 


[z-br»^^K{z,n). 


( 2 . 8 ) 


^ , , fe27vO°g(^ ^ +K(z,Q), if a is a right edge of an interval, 

4>(z) = 4>ao(-2) X < 1 , / N - 

|e 2 ,ri *°g(“ K{z, n), if a is a left edge of an interval. 

(2.9) 


Remark 2.5. The ‘h RH problem follows straightforwardly from the definition of 4>. 
However there is one slightly subtle point when computing the asymptotic behaviour 
in property (c). For in and z G C \ 71, where 71 is taken to be 71^®^®®) or ^(int) 
for edge and interior type problems respectively, the asymptotics claimed in (c) follows 
immediately. However if 2: —)• 00 for 2: G 71 \ 71;^ \ 71^ one still has the asymptotic 
behaviour in (jI.I2p or (jI.I3l) . with the associated K factor. However, the K factor can 
be removed by conjugating K through to the left and noting it contributes an error 
term smaller than 0(z~^). 


2.2 Examples 

In this section we will give some examples showing how the canonical model problem 
rednces to familiar model RH problems used previously in the random matrix theory 
literatnre. We make nse of the open lens version of the canonical model problem 
introduced above. 
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Figure 4: Jumps for the canonical model problem in the Bessel case at a hard edge 
with a logarithmic singularity at the edge. 

2.2.1 Relation to the Airy model problem 

Setting I = M, L = B = 0 and roo,o = 1 in a edge-type problem of order k = 0 we 
have the following RH problem, 

(a) $ : C \ M \ S® \ S® —)■ is analytic in z. 

(b) $ has the jumps shown in Figure [3] 

(c) As z ^ oo, 

$(^) = (/ + 0{z-^)) (2.10) 

This is exactly the Airy model problem |21j . 

2.2.2 Relation to the Bessel model problem 

Setting I = M“, L = M“, B = {0}, do = o/2 and roo,o = hr a edge-type problem of 
order A: = — 1, we have the following RH problem, 
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0 
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0 

-1 


1 

0 


0 
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Figure 5: The jumps for the canonical model problem in case of Bessel or Painleve II 
in the bulk with a logarithmic singularity. 


(b) 

(c) 

(d) 


<I> ; C \ M \ S* \ E® —)• is analytic in z. 

$ has the jumps shown in Figured 
As z ^ oo, 

$(z) = (/ + 0{z-^)) 

As z ^ 0, 

$(z) = 0(l)(z-6)t^3^(z). 

This can be written as, 


^ ^ 
Z 2 2; 2 


Hz) = 0 \ , 

V^;2 z 2 ' 

for 2 outside the lens and 


— d I _a _ a 

Vz 2 z 2 


for z inside the lens. 

This is exactly the Bessel model problem. 


( 2 . 11 ) 


( 2 . 12 ) 


(2.13) 


(2.14) 


2.2.3 Relation to Bessel in the bulk 

Setting / = M, L = M \ {0}, B = {0}, do = '^ 00,1 = 1 and roo,o = —^vra in an 

interior-type problem of order k = 0 we have the following RH problem, 

(a) : C \ M \ \ —>■ is analytic in z. 

(b) $ has the jumps shown in Figure [5] 

(c) As 2 ; ^ 00 , 

$(z) = (/ + Oiz-^)) (2.15) 
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The above RH problem is related to that in [22j by, 


■^{z) 


<h(2:) X 


^—■KioLB{z)az 

1 , 


if Re z > 0, 
if Re z < 0. 


(2.16) 


2.2.4 Relation to the general Painleve II eqnation 

Setting I = M, L = R \ {0}, B = {0}, do = a, Too ,3 = 4, Too ,2 = 0 and Too,i = 3s in an 
interior-type problem of order A: = 1 we have the following RH problem, 

(a) $ : C \ R \ \ —>■ is analytic in z. 

(b) $ has the jumps shown in Figure [5] 

(c) As z ^ 00 , 

$( 2 ) = (/ + 0{z-^)) e-3(4^"+3"^+^-.o)'"3Q(^), (2,17) 


The above RH problem is related to that in |19j by. 




(2.18) 


2.2.5 Relation to GUE at a birth of a cnt 

Setting / = R, R = 0, Too ,2 = Ij Too,i = 0 and Too,o = 0 in an exterior-type problem of 
order A; = 1 we have the following RH problem, 

(a) $ : C \ R —>■ is analytic in 2 :. 

(b) <i> has the jump, 

$-H2)$(2) = I) (2.19) 

for 2 : G R. 

(c) As 2 ^ 00 , 

$(z) = (7 + 0(2"^)) (2.20) 

This RH problem is related to the standard RH problem for Hermite polynomials 
of order c by, 

$( 2 ) = y( 2 )e"^'"T (2,21) 


2.3 Solvability of <|) RH problem: Proof of theorem 11.241 


We follow the standard argument used in [23] . In particular the existence of a solution 
to the model Riemann Hilbert problem depends crucially on showing a “vanishing 
lemma”. To state it, we define a function A so that we may write the asymptotics 
()1.12l) and (I1.13P in a uniform manner. In particular let. 


\-^,p2k+i{z)e{z) 

we then have the following result. 


if is of edge type, 
if <I> is of interior type. 


( 2 . 22 ) 
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Lemma 2.6 (Vanishing Lemma). Let $0 be a funetion satisfying the conditions (a), 
(b) and (d) of the RH problem for $ together with the asymptotic behaviour as z ^ 00 
of, 


$o(^) = n)-^ (2.23) 

where O takes the form specified in the RH problem for in the edge and interior eases, 
r] = 2,!A if ^ is of edge type and r] = 1 if ^ is of interior type. Then $0 = 0. 

Proof. Let us first draw the readers attention to the following properties of A. Let A{x) 
be a function which is real for a: G M. When <1> is edge type we have, 


A+(a:) = A_(x) = —A(x), x > 0 , (2.24) 

A_|_(x) = —A_(x) = —fA(x), X < 0, (2.25) 

and when < 1 > is of interior type we have 

A_|_(x) = —A_(x) = —fA(x), X G M. (2.26) 

In the above we have used the fact that the vector Too is real. 

We now let, 

G{z) := $o(^)e-^(")'^®Q(z) (2.27) 

where we have used, 


Q{z) 



Im z > 0 , 


I Imz < 0. 

The function G satisfies the RH problem. 


(2.28) 


RH problem for G 

(a) G : C \ M \ —>■ is analytic. 

(b) The jump matrix jciz) '.= G-{z)~^G^{z) has the following form, 

Xl{z)e^+ R) —e^- (z) 

g-A_(2:)+A+(2)-27riarH) Q 



if z G 


(2.29) 

jc = e^ 2 «ar( 2 )a 3 ^ if > 0 and G Uftes (2.30) 

These jumps specialise in the edge and interior cases to: 

(i) < 1 > is of edge type: 


JG 


JG 


g —g27r2Qr(2) 

g-27rjarR) Q 


z G M+ 


Xl{z) _^2iA{z)+2m&Y{z) 

^ — 2iK{z)—2'Kia.Y{z) g 


z G M 


(2.31) 

(2.32) 


JQ = eG27riar(H<^3^ ff ^ > 0 and z G Ubes Tfo. (2.33) 
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(ii) $ is of interior type: 


JG 


e 


Xi{z) 

2 iA{z)— 27 riaY {z) 


^2iA{z)^2TTia-piz) 


0 


z G M 


JG = if > 0 and z G IJfeeB 


(2.34) 

(2.35) 


(c) As z ^ oo, 

G{z) = 0{z-^), (2.36) 

where r/ = 3/4 if ‘h is of edge type and r; = 1 if 4> is of interior type. Note that 
this behaviour depends crucially on the constraint roo, 2 fc+i > 0 in the interior 
case, roo,fc+i > 0 in the edge case with k G 2N*^ and Too,fc+i < 0 in the edge case 
with k = —1. 


(d) As z —)• 6 G B, 


G{z) = 0(1) exp 


di,-l 


1 


i=i 


1 log \ z — a\ 

1 log \z — a\) ’ 


As z ^ a £ dl \ B for Im z < 0 we have, 
G{z} = O 

while for z ^ a £ dI\B for Im z > 0, 
G{z} = O 


log \z — a\ 1 
log \z — a\ 1J ' 


[z - bf^'^^Qiz) 


We now introduce the function, 

H{z) :=G{z)Giz*)\ z G C \ M. 


(2.37) 


(2.38) 


(2.39) 


(2.40) 


Using the jumps of G it is indeed possible to demonstrate that H has no jump on 
Ubg^rf, \ M and no singularities at b £ B and is therefore analytic in C \ M. The only 
points at which some elements of H are not bounded are in dI\B, at which there may 
be log divergences. Lastly, we have that as z —>• oo, 

H{z) = 0{z-^^). (2.41) 


Using the above properties together with Cauchy’s theorem allows us to conclude, 


H^{x)dx = 0. 


(2.42) 


By adding the hermitian conjugate of the above equation to itself we have, 


L 


/n(—oo,o] 


G-W(J “iG- 


(xydx+ f G_(a:) O G_(a:)^da: - 0, (2.43} 

4/n[0,oo) \ U uy 
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where V(x) := in the edge case and V(x) := 1 in the interior case. In both cases 

we are able to conclude that the first column and second column of G(z} is zero in the 
lower and upper half plane respectively. 

To analyse the other entries of G we use the standard argument based on Carlson’s 
theorem [23| . We dehne the scalar functions, 


Gk,i(z) lmz>0 
Gk, 2 (z) Im2;<0 


(2.44) 


and it is then easy to show that both gi and <72 satisfy the following RH problem: 


RH problem for g 

(a) g : C\R Uh^B F/, —>■ C is analytic. 

(b) g has jumps, 

g+iz) = ^_(^)e- 2 *AC)Vi(^)- 2 ™r( 2 )^ ^ (2.45) 

g+{z) = 5 -(^)e-2™rC), z G U TA K (2.46) 

beB 

where Vi(^) = Xk-('2 ) in the edge case and Vi{z) = 1 in the interior case. 

(c) As 2 : —)■ 00 , g{z) = 0{z~'^). 

(d) As 2 —)• 6 G R we have. 


g{z) = 0(1) exp 
As z ^ a G (9/ we have g{z) = 0(1). 


du — l 


i=i ■' 


[z-b] 


-«6 


(2.47) 


Remark 2.7. In condition (d) we have g{z) = 0(1) as 2 ; ^ a G dl. In principle we 
could have had g{z) = 0{log{z — a)), however it is easy to see by deleting the jumps 
near a that we must have 0(1) behaviour at a. 

We now transform g in order to put it into a known form. Define, 

g{z) := g{z)z~^^°^ 11 ( ~ 

h&B \ 

We then have, 






i=i 


(2.48) 


RH problem for g 

(a) ^ : C \ M is analytic. 

(b) g has jumps, 

Y{z) = 2 G M (2.49) 
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(c) As z ^ oo, g{z) = 0{z 

(d) As z —0, g(z) = 


At this point we note that for the edge case a problem of this form has been solved in 
|23j and for the interior case this form of problem has been solved in m- In both cases 
one finds, via Carlson’s theorem, that g = 0. The only difference in the current case 
compared to [23l[T9] is that here A is not explicit. Nevertheless, the form of proof used 
in m does not rely on the precise form of A and therefore the result follows immediately 
in the interior case. For the edge case we have essentially the same situation, however 
we feel it is useful to highlight some of the changes necessary in the proof from |23j . 
We first follow [23] by making a change of variable. 


h{z) 


g{z‘^), Re z > 0 

, Re ^ < 0. 


(2.50) 


It is straight-forward to verify that h is analytic in the region Re 2 ; > 0. The difference 
compared to [23j arises in the next transformation. 


h{z) 


1 + z 


4k±8. 

2k+3 


2fc+4 
2; 2fe+3 


(2.51) 


Note that the above equation reduces to equation (2.29) in [23| in the case k = 0. One 
can then verify that h is analytic for Re z > 0, is bounded for Re z > 0 and for z —)■ 00 
on the line Re 2 ; = 0 we have, 

\h{ix)\ < (2.52) 

where C,c > 0 are constants. By Carlson’s theorem we are able to conclude that h = 0. 
This completes the proof of the vanishing lemma. □ 


The proof of Theorem 11.241 follows from the description of the RH problem in 
terms of singular integral equations of Cauchy-type whose corresponding operator is 
a Fredholm operator of index zero. It can be shown that the kernel of this operator 
is trivial if and only if the vanishing lemma holds and therefore the vanishing lemma 
implies the integral equation is solvable. A detailed description of these points can be 
found in |23| . One subtlety of this argument is that the RH problem under consideration 
must be equivalent to one with no singular points which is not the case here. On this 
point we follow the same argument given in the existence proof found in |15j . 

Noting the asymptotic behaviour of ^ as 2 —)• 6 € R, we make the following trans¬ 
formation <I> I—)• <I>fe o(2:), 


^{z) = 4>b,o(2;)exp 


db-l 

E 

j=i 


^Tlj{z-b) V 3 


[z-b] 


060-3 


(2.53) 


where <I';,^o(-2) is bounded as z ^ b £ B. Substituting the above relation into the jump 
conditions yield that has no jumps in a neighbourhood of 6 if 6 ^ M. This implies 
that is analytic near such a point. For 6 € M we instead find that must satisfy 
the jump jo := ^b,o-{z)~^^b,o,+ {z), 


jo 


g2vri(Qr(2)-XrJ^)o6) 

0 


XI (z) 


_ ^ 1 . ^ 
p 7 = 1 70 . 


.1 ■ (z—b) ^ I 71 

J = J I I 2 : _ 5| 


3 0,3 

-2iTi{ar (z)-Xr^ (^job) 


2o6 


ZGIU [J Ffe/. 
b'&B 

(2.54) 
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We now let, 


(2.55) 






Rz) 

1 


g7ri(or iz)-Xr^ (z)6n,)d(z)a3 


and note that if we choose / such that, 


f+{z) = f-{z)+xi{z)e ^^=1 


\z-b\ 


2o(, 


(2.56) 


for z in any fixed interval around b we have that ‘hft q has no jumps in a fixed neigh¬ 
bourhood of b and is hence analytic there. It is important to note at this point that a 

V ^ ^ 1 j / 7 \ _ j ^ 

solution of (j2.56p only exists if is such that e~> |^_ 5 | 2 a 6 jg bounded 

on I. This leads directly to the constraints (|1.53p stated in the theorem. 

Let Dh be a fixed disc centered on b G B. We now define. 


^(z) = 


$(z). 


‘^6,o(^) UhPR^b 


z gC \ Ufees ^b, 
b&B ^ 


(2.57) 


and note that has no singular points. This completes the proof for the edge and 
interior case. In the exterior case we postpone a proof until Remark 13.51 

Remark 2.8. We have shown that the canonical RH problem is solvable for a certain 
set of admissible real vectors Too- The transformation to a RH problem with no singular 
points also lets us use the same argument as used in [23] to show that the canonical RH 
problem is solvable for complex parameters in a neighbourhood of such an admissible 
vector. 


3 Asymptotic analysis of the RH problem for orthogonal 
polynomials 

3.1 The RH problem for orthogonal polynomials 

An effective way to characterise orthogonal polynomials appearing (11.41) is via a well 
known RH problem due to Fokas-Its-Kitaev |25j . 


RH problem for Y 

(a) Y : C \ X —^ is analytic. 

(b) The limits of T as z approaches M from above and below exist, are continuous on 
M and are denoted by 1+ and IT respectively. Furthermore they are related by 

y+(a;) = y_(a;)(j , x gX. ( 3 . 1 ) 


(c) Y(z) = {I + 0{z ^))z”'^3 as z ^ oo. 

(d) As z —)• a for o G .4. \ B, 


Y(z) = Ya{z) X 




if a is a right edge of an interval, 
if a is a left edge of an interval. 


(3.2) 
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where is analytic in a neighbourhood of a. As z ^ b for b € B we have, 


Y{z) = O (I 



(3.3) 


Remark 3.1. The behaviour (13.21) follows from the fact that Y has at most a log 
divergence at a and Y has the same jumps as exp(^'u;(z) log(z — a)cr+) at a right 
edge. A similar statement holds at a left edge. 

This RH problem has a unique solution, 


y(z) 


f Pniz) qn{z) 


where pj is the degree j monic orthogonal polynomial defined in (II.3p and 


Qjiz) 


J_ r pj{x)w{x) 
2TTi Jx X — z 


(3.4) 


(3.5) 


3.2 Transformation to constant jumps 

To construct the first transformation we introduce an analytic continuation of tCbr- 

Definition 3.2. We define. 


Wh,{z) := 

beB 

Let w{z) := Wreg{z)Wsing{z)Wbriz). 

Remark 3.3. We have that tcbr has the following jump properties. 


Whr{z)+ = Wbr( 2 :)-e"‘™r^ 


(3.7) 


beB 


Whriz) + Wbriz)- = 

z gR. 

(3.8) 


fees 


Note that the second jump property relies on the fact that B contains conjugate pairs 
of points. 

Definition 3.4. Define or(- 2 ) in an identical way to Definition 11.71 with B replacing B 
and Ofo replacing ai,. 

(73 

Defining 'I'(z) := Y{z)w{z) 2 , we have. 


RH Problem for T 

(a) T : C \ X \ UfeeB ^6 —>■ is analytic. 

(b) Let j^{z) := 'I'_(z)“^'I'+(z). Then, 

J'^yZ) I Q g-2i7ror(2) J ’ 


zeXuljTfe. (3.9) 

beB 
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(c) '£'(- 2 ;) = (/ + 0(z ^ as z ^ oo, where atot := YlbeB^b- 

(d) As z —)• 6 s B, 


'I'( 2 ) = 0(1) exp 


db-l 


iZ ^^b,jiz-b) Vs 


i=i 


[z-b] 


ctbcrs 


(3.10) 


As z ^ a £ A \ B 




^a,o(^) X 


e-^i^°s{a.-z)a+ ^ 


if a is a right edge of an interval, 
if a is a left edge of an interval. 


(3.11) 


Here 'ha,o is analytic with respect to 2 : at a. 


Remark 3.5. Observe that the RH problem for exactly matches the definition of 
the canonical model problem of exterior type and we therefore are able to construct a 
solution to a model problem in terms of orthogonal polynomials. 


3.3 Opening the lens 

We now perform the standard step of opening the lens. The difference here is that, 
since the jumps on the lens contours are constant, the lens contours are unconstrained. 
In a subsequent transformation we will use the g-function to normalise at infinity, at 
which point the lens contours will be required to stay within a region in which they 
converge to the identity as n ^ 00 . 

Due to the rather general nature of the problem we need some additional definitions 
in order to define the lens contours and associated regions. 


3.3.1 Definitions of contours 

Definition 3.6. For an interval cr of M define the contours to be smooth contours 
from inf a to sup a in the regions a ± zM'*' respectively. A lens contour for a is the 
contour U T,Z_ with all contours orientated from inf a to sup a . 

Definition 3.7. Given a lens contour for an interval a define to be the region of 
C bounded by cr U such that the + side and — side of a is contained in and VLZ. 
respectively. 

Definition 3.8. For a collection of pairwise disjoint intervals X define a full lens contour 
TA as all intervals in X treated as contours together with their lens contours. Denote 
the lens contours of YA by and define := IJo-ex ^±' 

Finally we surround all singular points with discs in which we will later construct 
the local parametricies. 

Definition 3.9. For all x* G P* U 7^* define, Dx^(S) to be a disc of diameter delta 
centered at x* whose boundary is orientated clockwise. Let U = Ux*e 7 ^*u 7 e* 

We will always work with n large enough that the discs contain all points in P U 7?.. 

Let S C S nU he the intervals of S which are fully contained in U, i.e. they don’t 
intersect dU. Let 7 = 5 \ 5 \ 7^*. We now open the lens of 7 by defining, 

5 ( 2 ) := ^(z)K(z), ( 3 . 12 ) 
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Figure 6: An example of contours constructed using Definitions 13.6113.71 and 13.81 


where AT is a piecewise function designed to open the lens, 


(I, ior z £C\nl\Ql, 

K{z) := \^_e2™r(D i) ’ for ^ e 

g-2iar(D J) ’ for z G DI. 

The function S satishes the following RH problem. 


(3.13) 


RH Problem for S 

(a) 5 : C \ X \ ST" \ (JfeeB F^ ^ is analytic. 

(b) Let js{z) := S'_(2;)“^5+(2;) then, 


Js = 


JS = 


0 1 
-1 0 y ’ 

g2ivrar(D 

^—2i7Ta-p{z) I ’ 


_ 1 0 

Js — g±2i7rar(2) J I ’ 


2 G 7, 

z GXU IJ Fb\7, 

beB 


z G SX- 


(c) S{z) has the same asymptotic behaviour as '^(z) as z ^ oo. 


(3.14) 

(3.15) 

(3.16) 


(d) S{z) has the same asymptotic behaviour as ^'(z) as z ^ 6 for 6 G .A U R unless 
we also have 6 G 5, in which case. 


S{z) = ^{z) X < 


for z G C \ \ 


1 

o\ 


g27ria:r(2) 


, for z G 

1 

o\ 

for 2 ; G 

—27rior ( 2 ) 

lb 


(3.17) 
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3.4 Normalisation at infinity 

The next transformation takes the form. 


^ ifzec\u, 

T{z) := e 2 S{z) x 

I, if z G 17. 


(3.18) 


The above transformation has the effect of normalising the problem at infinity. 


RH Problem for T 

(a) T : C \ X \ XT' \ (91/ \ Ufeeg ^ ^ analytic. 

(b) Let jT(^) ■= T-{z)~^T^{z). Then, 

JT = js, 


JT 

JT 


0 

-1 


1 

0 


^2Txi{aY{z)-ne{z)) 

0 


^-2-Ki(aY{z)-ne{z)) 


JT = 


1 

^■^2'KioLx^z—2n^{z) 


.JT = 


(c) As z —>• oo, 

T{z) = (l + 0(z-^))z"*°‘. 


zG (x\s^\ |Jrfe)nt/, 

beB 

(3.19) 

ze-f\U, (3.20) 


2; G X U 

n\i\u, 

beB 

(3.21) 

zgj:1\u 

(3.22) 

z G dU. 

(3.23) 


(3.24) 


(d) As z ^ b G AU B, T{z) has the same asymptotics as e 
for sufficiently large n all points in .4 U 13 are in U. 


nla^ 

2 S{z). This is because 


3.5 Final transformation 


For the final transformation we need a global approximation to the RH problem for T. 
To define this we first need to introduce some new objects. 


Definition 3.10. Define := e(x* ± ^^x*) and := ar(3i* 
the values of e and or at the points at which dDx^, intersects 


± 


these are 


We also define 


;= £ — 

X* X:ii 

and ai 

:= — at ■ Note that et 

can be undefined when x* is of 

edge or interior type. 



Definition 3.11. Let, 




ar(2:), 

for 2 ; G M \ 17, 


Q;r(^) = < 

a"*' 

for 2; G Dx^ n (x*, 00), 

(3.25) 


a~ , 

V. X* 5 

for 2; G Dx, n (—oo,x*). 
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Figure 7: On the left the figure shows the behaviour of e on the connected components 
of M \ (shown as dashed lines) in the case that x* is an exterior point. The value of 
e at the boundary of the disc are labelled In the figure on the right we show the 
behaviour of ei introduced in Definition 13.121 in which the values of e on the boundary 
of the disc are extended all the way to x*. 


Definition 3.12. Let, 


h{z) 


e{z}, for z G M \ 5* \ t/, 

< for 2 G Dx^ n (x*,oo) with x* an edge or exterior point, (3.26) 

e“^, for z G Dx^ H (—oo,x*) with x* an exterior point. 


Note that ii is constant on connected components of Now we define e{z) := 

ei{z) — n~^M{z) where : M \ JL, —)• N*’ is a function which is constant on each 
connected component of M \ ^7*. We chose M{z) such that when x* is an exterior point 
we have n\ix^ “ ^ = 0 if nei{z) G Z. 

Definition 3.13. Define := e(x* ± ^d^*) and := 


Remark 3.14. The idea in the definition of dr and e is that we want them to match 
ar and e where they are defined but extend them to the domain obtained by shrinking 
the discs in U to points (see Figure El). The reason for this is that we can then use 
a simpler global parametrix built on a Riemann surface which is not degenerate as 
n —oo. Furthermore, in the definition of e we have used the freedom that the jumps 
of T are invariant under ne^ ne + k where k is integer, to minimise the difference in 
e across singular points x*. 

Finally we define the global parametrix as a function satisfying the following 

RH problem, 

RH problem for 

(a) : C \ (—oo,sup(.S* U J^)] is analytic. 

(b) has the jump relations 

p|-)(^) = pM(^)(^_0^ q), zG5, (3.27) 

{z) = (z)e2^*(“r (D-ne(D)^3 ^ z G M \ J*, (3.28) 

Recall that dr : M \ —)• M and e : M \ ^7* —t M are functions which are constant 
on each connected component of R\i7*. Furthermore e is zero on infinite intervals 
and dr is zero on M \ (—oo, sup.8*). 
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(c) As z ^ oo, 


p(^\z) = (/ + 0(2-^))2“*°‘^L 

(3.29) 

(d) As 2 ^ 2 G 55* 


pM(^) = C)((z-.2)-i). 

(3.30) 

As z —)■ G 


pi°^){z) = 0{l)iz - 

(3.31) 


We now need to define a local parametrix for x* G P* U 7^*. We define 

as the solution to the following RH problem, 


RH problem for 

(a) : R>a:* (4.) \ 2 : \ XT' \ UfeeB^b ^ is analytic. 

(b) Let jp(:n,) := P^*\z) then jp(.*) = jt- 

(c) As z ^ b G B\x,, P^^*\z) = T{z) where the equality is understood to mean equal¬ 
ity between asymptotics series and the notation B\x, was introduced in Definition 

(d) As n —)• 00 , P^^*\z) = {I + uniformly for G dDx,{Sx,)- 


We now define R{z) as, 


R{z) 


r(^)p(oo)(^)-i^ 

T{z)P^^*\z)-\ 


for z G C \ [/, 

for 2 : G Dx,{6x,), where x* G P* U 7^*. 


(3.32) 


Using the above definition together with Definition 13.121 we can derive the following 
RH problem for R. 


RH problem for R 

(a) R:C\I\iJ:'^\U)\dU ^ is analytic. 

(b) Let jniz) := R-{z)~^Rp{z) then, 

jn(z) = Pl“)(z) (^_ 2 .^,,,L™p(.) () 

(\ pfi(?-(2)+S+U))-l-27rme(z)-27riar(2)\ 

iHW = pi“’w(„ j jpUw-', 

■ ( \ x>P*)( 0 \ (oo)p x-i 

3 r{z) = P_ {z) ( ^ 1 P+ {z) , 


(c) As 2 —)■ 00 , 7 ^( 2 ;) = / + 0{z ^). 


2 G SX \ U, 

(3.33) 

z€l\j\U, 

(3.34) 

2: G dDx, , 

(3.35) 
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3.6 Proof that R is a small norm RH problem 

To prove that R is a small norm RH problem we will show that each jump matrix of R 
tends towards the identity as n —)■ oo. The first observation is that for sufficiently large 
n and any fixed point x & S \ U we have that p{x) > 0 |21j . By standard arguments 
this implies that there exists a choice of lens contours such that Re.^(z) > 0 for 
all z G \ U which together with the form of the jump matrices implies the jump 
matrices converge uniformly to the identity on \ U. 

We now turn our attention to ()3.34l) . In this case we make use of the variational 
condition that ^+(x) + ^-(x) < 0 for z G M \ 5 \ C/ to again conclude that the jump 
matrices converge uniformly to the identity on M\S \ U. 

Finally we turn our attention to the jump on the disc boundary (I3.35h . That this 
jump tends to the identity as n —>• oo follows trivially from property (d) of the local 
parametrix. 

We therefore conclude that all the jumps for R decay to the identity as n —)• oo. 


4 Construction of global and local parametrices 

4.1 Construction of 

Given a collection of closed pairwise disjoint intervals I we hrst define a scalar Szego 
function D{z;a,I) that solves the following RH problem, 


RH problem for D 

(a) D{z; a, I) is non-zero and analytic on C \ M with respect to 2 ;. 

(b) D{z;a,I) satishes the following jump relations: 

T)+(x; a,/)il_ (x; a,/) = for x G/, (4.1) 

T)+(x;q;,/) = for x G M \/, (4.2) 

where the function : M \ / —)• M is constant on each connected component of 

M\/. 

(c) D and D~^ remain bounded as z ^ z £ dl and 

Dao ■= lim D{z-,a,I) (4.3) 

z—>-00 

exists and is non-zero. 

(d) As 2 ; —0 we have, 


D{z) 


z"(l + C>(z)), ifOGl 
0(1), ifO^/' 


(4.4) 


Remark 4.1. The function ija in the definition of D is not arbitrary and is fully 
determined by the requirement (c) in the RH problem. It is shown in |22j how (c) 
determines r]a explicitly, together with an explicit construction of D. 
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Definition 4.2. We now define two auxiliary functions built using D; 


(4.5) 

(4.6) 


D{z;x^:,a) := 
D{z-,x^,,a) := 


(z —X*) "^(z — X*; a, 5* 
D{z]bo,a) 


D{z; X*, a) 


X* ^ iS*, 


X*) 


where 6 o := sup(5* H (—oo,x*)) and principal branches are taken for multi-valued 
functions. The notation 5* — x* means the intervals of 5* shifted to the right by x*. 

It is straightforward to verify that D and D satisfy the following RH problems, 


RH problem for D 

(a) The function l) : C \ M —)■ C is analytic in z. 

(b) D has jumps, 


b+{z)b.{z) = 1 

z G 5* 

(4.7) 

b+{z) = D_(z)e2"^(^), 

z G (x*, 00) n (M \ 5 *), 

(4.8) 

D+(z) = D_(z)e2"(^(^)-0, 

2: G (—00, X*) n (M \ 5 *), 

(4.9) 


where 77 : M \ 5* —)• K is a function constant on each connected component of 
M \ 5* and zero on infinite intervals. 

(c) As z ^ 00 , 

D{z) = z-^D^{l + o{l)). (4.10) 


(d) As z ^ X*, 


D{z) 


{z — x^) “ 0 ( 1 ), if X* ^ 5* 
0(1) if X* G 5* 


(4.11) 


RH problem for D 

(a) The function D : C \ M —)• C is analytic. 

(b) D has jumps. 


b^{z)D_{z) = l 

2 G 5* 

(4.12) 

b+{z) = 

2 G M \ 5* \ [ 60 ) 

(4.13) 

b+{z) = 

2 : G [& 0 ) a:*] n (M \ 5*), 

(4.14) 


where 77 : M \ 5* —^ M is a function constant on each connected component of 
M \ 5* and zero on infinite intervals. 

(c) As z ^ 00 , 

D{z) = D^{l + 0{z-^)). (4.15) 
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(d) As z ^ X*, 


D{z) = (z-x*)“0(l). 

(4.16) 

As z ^ bo, 


D{z) = 0{l). 

(4.17) 


Note that the above behaviour is obtained under the constraint appearing in the 
definition of D that x* ^ 5*. 

We now seek in the form, 

n J] D^{b, x 

x*e£* feeB* 

n D{z-,b,abr n D{z-x.,neij\ (4.18) 

Note that all points in and are real. Furthermore, using the definition of ar{z) 
(Definition 13.111) together with the fact that e — is constant on 

connected components of R \ ^7* we obtain that the function must solve a RH 

problem of the form, 

RH problem for 

(a) : C \ R —)■ is analytic in z. 

(b) has the following jump relations: 

P^\x) = P^\x) , for a; G 5*, (4.19) 

pM(x) = p!°°)(x)(^ Q forxGR\5,. (4.20) 

where v{z) : R \ 5* —)■ R is a function constant on each connected component of 
5* \ R and zero on infinite intervals. 

(c) As z —)■ oo, P^°°\z) = I + 0{l/z). 

The solution to the RH problem for in terms of theta functions is well known and 
can be found in |22j . 

4.2 Asymptotic behaviour of ^ and near singular points 

Before constructing the local parametrix we first prove some useful lemmas regarding 
the behaviour of and ^ near x*. We will make heavy use of these lemmas in the 
construction of In the following we consider the case of x* being an interior, 

exterior or right edge point. The case of a left edge point we omit as it is similar to the 
right edge case. 

Lemma 4.3. Consider P^°°\z) inside the disc Dx,- It will be convenient to change to 
the variable z = x* + C G dDx^ ■ Let P(C) be an analytic function in a fixed neighbour¬ 
hood of C = 0- 
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(4.21) 


(i) If G 5*, i.e. an interior point, then we have, 

p(-)(x, + c) = ^(C)Q(C), 

where we have used Q as defined in (I1.14p . 

(a) If G dS^ with X* forming the right edge of an interval in 95*, i.e. a right edge 
point, then we have, 

P^°°Hx^ + C) = (4.22) 

where we recall N := 2 ~2 (^I -\- iai). 

(Hi) If x^ £ 5* then we have, 

P(°°)(x* + C) = (4.23) 

where we have introduced (3 := — ni^^. 

Proof. We consider each case in turn. 

(i) If X* G 5* then E{z) := P^°°\z)Q{z)~^ has no jump in a neighbourhood of x*. 
Furthermore E{z) is bounded at x* and therefore we conclude E is an analytic 
function in this neighbourhood. Redefining E by shifting the argument we obtain 
the result (i). 

(ii) If X* G 95* with X* forming the right edge of an interval in 95*, then, 

E{z) := _ x*)^ (4.24) 

has no jumps in a neighbourhood M of x* and therefore E is analytic in M\{x*}. 
Furthermore, given that the behaviour of P^°°\z) is 0((z — x*)“4) near x* we 
see that E can be analytically continued to all of M. Redefining E by shifting 
the argument we obtain the result (ii). 

(iii) If X* G 5* then, 

E{z) := (4.25) 

where /3 is defined in the statement of the theorem, has no jumps in a fixed 
neighbourhood of x*. Given the behaviour of P^°°\z) as 2 ; —>■ x*, ()3.3ip . we see 
that E is analytic at x*. Redefining E by shifting the argument we obtain the 
result (iii). 

□ 

Remark 4.4. For a left edge one finds, 

F(°°)(x* +C) = (4.26) 

Lemma 4.5. Consider (^(z) on the boundary of a disc Dx,. The points on the boundary 
of the disc may be parameterised by letting z = x* + (" with |C| fixed. Let Ei for I G Z 
be functions analytic in a neighbourhood of ( = 0 which behave as 0{1) as n —)• 00 
uniformly for Q in a neighbourhood of zero. We then have, 
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(i) If G 5* is an interior point of order k, then as n —)• oo 


+ C) = <+(a:*)0(C) - J](n^-C)''=+^"'^2A:+i-KC)0(C), (4.27) 

where, for sufficiently large n, £' 2^+1 (0) > 0 and E 2 k+i-i{ 0 ) are real. 

(a) If G 95* is an edge point of order k, with x* forming the right edge of an 
interval in 95*, then as n —)■ 00 , 


<(x* + C) = 7rm6l(C)e+ 

where, for sufficiently large n, 
and Ek+i{0) <0 if k = —1. 


00 

1=0 

Ek+i-i{0) are real and £fc+i(0) >0, if k ^ 2N° 


(in) If x^: G 5* is an exterior point of order k, then as n ^ 00 , 


n^(x* + C) 


7rm6l(C)e+ +ne^^ log (n^"*C) 


1 


/=o 


(4.29) 


where for sufficiently large n, E 2 k-i{ 0 ) are real. 

Proof. The proof will proceed by establishing the existence of a Laurent series represen¬ 
tation for f on dDx^. We will do this in two ways; firstly using the jump properties of 
f and then by using (jl..S5p . Both methods give us some, but not all, of the information 
required and we obtain the result by combining these properties using the uniqueness 
of the Laurent series. We now establish the basic form that the asymptotics of f must 
take in each of the following cases. 


(i) X* is an interior point. Define f,{z) := f^{z)6{z) for z G D^*- Note that for 
sufficiently large n, ^ has no jumps that intersect the boundary of the disc which 
together with Assumption 11.191 implies there exists some fixed annulus centered 
at X* in which ^ is analytic. Choosing dDx, to lie within this annulus we have 
that on dDx„, ^ can be written as a convergent Laurent expansion centered at 
X*. This gives, 

OO 

^(x* +0 = 0(0 fjC (4.30) 

j=-ao 

where are Oindependent constants. We have yet to specify the asymptotic 
behaviour of 0 as n —00 but (I4.30p implies that the large n asymptotics of 
0^* + 0 must have this functional dependence on 0 

(ii) X* is a right edge point. Define f,{z) := f{z) — iTief^9{z)-, this removes one of the 
jumps that extends to dDx„. The other jump extending to dDx, can be removed 
by defining 0 (. 2 ) := f,{z){z — x*)“^/^ where the principal branch is taken with a 
cut on M“. Again, because dDx^ lies in the annulus in which 0 is analytic we 
conclude that. 


OO 

O3^*+o = 7r+O0(o + c^ CjC^- 

j=-oo 


(4.31) 
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(iii) X* is an exterior point. Define .^(z) := ^{z) — TTie^^6{z); this removes one of the 
jumps that extends to dDx^- The other jump extending to dDx, can be removed 
by dehning ^i(z) := ^{z) — (e“^ — log{z — x*). As before we can then write, 

OO 

C(x* + C) = vrfe+ 0(C) + {e-^ - e+J logC + ^ (4.32) 

j = -QO 


We note at this point that we are able to rearrange M.SOp . (I4.31h and (I4.32h into 
the form ()4.27l) . ()4.28l) and ()4.29p respectively. What is missing is the n dependence 
of the coefficients in the Laurent expansion and the properties of Ei. To obtain these 
details we now use (ll.35p . First we write it as, 

?(^) = i±{v') j ^ y{s)ds, ±lmz > 0, (4.33) 

where the integration contour does not intersect M and p' = sup(77.|a;*), i.e. the right 
most interval edge in the disc, if sup(77.|a:*) exists, otherwise we set p' = x*. Note that 
when X* is of right edge or exterior type we have ^±{p') = so we may write, 

^{z) = 7rf0(z)e+ - y{s)ds, (4.34) 

whereas when x* is an interior point we have, ^±{p') = ±^+(4*0 so we may write, 

C(2) = 0(z)C+(p') y{s)ds. (4.35) 

Now note that y has a jump relation, 

y+{x) + y-(x) = 0, X G 5. (4.36) 


By mimicking the arguments used previously for C we may delete the jumps inside 
and write. 


y{z) = hi{z)y^Ri{z)E{z) x 



X* is a exterior or edge point, 
X* is an interior point. 


(4.37) 


Here hi is a polynomial and Ri is a rational function with only simple poles, with both 
functions only having poles and roots in Dx^- The function E[z) is analytic in a fixed 
neighbourhood of z = x* with no zeros in Dx^ ■ 

Remark 4.6. We note that the factors hi, Ri and E are only defined up to overall 
constants and we have some freedom how we normalise each factor. We choose E[x^) = 

1 . 

In light of Assumption 11.191 we have that if hi and Ri have ruh and rriR zeros 
respectively and Ri has nip poles we can write, 

/ii(x* /i((t) (4.38) 

Ri (x* + a) = R{a) (4.39) 
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where R{z) = ri{z)/r 2 {z) and fi, r 2 and h are polynomials in which each root behaves 
as 0(1) as n —)• oo. Note that since the behaviour of lim„^ooy('^) near x* takes the 
form in Remark [1.131 we must have (nih + + 1)A2;^ = 1. 

This all leads us to define, 


yi{z) := h{z)\jR{z), 
which we note is related to (|1.45l) by, 


lim x 

n—)-oo 

= lim yi(C) = yiO- 

n^oo 



X* is a exterior or edge point, 
X* is an interior point. 


(4.40) 


(4.41) 


Define p = (p' — x*) and note that this also behaves asO(l)asn—>-oo. Using 

the taylor series for E(z) about x* we then have. 


-j OO 

C(x*+C) = C±(p )- 9 X] 

j=o 


C, 


yi(a)a^ dax 


1 , 

0(0, 


X* is a exterior or edge point, 
X* is an interior point. 

(4.42) 


where the coefficients ej of the taylor series of E have 0(1) behaviour as n —?> oo. 

To compute large n asymptotics of the expression (I4.42j) we need large ( asymptotics 
for the function, 

yi((T)cr^da. (4.43) 

Such asymptotics for Yj can be obtained simply by integrating the large C behaviour of 
yi term-by-term and adding an arbitrary constant; this is valid because yi(C) can be 
written as a convergent series for large C- 

We now again consider each type of critical point in turn. 

(i) When x* is an interior point, we have as C —^ oo, 

OO 

= (4.44) 

1=0 

for /c G N^, where qi are constants. We therefore have as C, oo, 

OO 

r,(c) = -^^c“+‘+^Ewr', (4.45) 

1=0 

where qj^i := (2k + l)/(2i(2k + j — 1 + l))qi are constants. Note that the above 
expression could have in principle contained a logarithm, however this would have 
been inconsistent with (j4.30p . Using the above expression we obtain. 


<(x* -h C) = 
<+(x*)6l(C) 


2k + 


y0(C) 


/=o 


3=0 


(4.46) 


where we have changed ^+(p') to ^+(x*) by shifting the value of q 2 k+i,o- Define 
E2k+i-i(C) ■= Yl^o^ 3 'l 3 ,lC^ ^*4 note that, using the expression for qj^i, it is 
absolutely convergent in a neighbourhood of zero and therefore analytic there. 
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(ii) When x* is a right edge point, we have as C —^ oo, 

OO 

yi(C) = (4.47) 

1=0 

where k G {—1} U 2N®. We therefore have as C oO; 

oo 

yjic) = ( 4 . 48 ) 

1=0 

where qj^i := {2k + 3)/{2{2k + 2j — 2/ + 3))^; are constants. Note that the above 
expression could have in principle contained a constant term, however this would 
have been inconsistent with (I4.3ip . Using the above expression we obtain. 


+ C) = 


Define Ek+i-i{C) ■= 
function at zero. 


(4.49) 


2 

2A; + 3 


1=0 


j=0 


E OO 

j=0 


ejqj,iC^ 


and note that again it defines an analytic 


(iii) When x* is an exterior point, we have as C —>• co, 

CX) 

MC) = 

1=0 

for A: G N. We therefore have as C ^ co 

1 1 °° 

U(C) = -:^»+jlog(0 - «C“+'E9 jX“'. 


(4.50) 


(4.51) 


1=0 


where qj^i = 2k/{2k + j — l)qi for I ^ 2k + j and qj, 2 k+j are constants. Using the 
above expression we obtain. 


n^(x* + C)= (4.52) 

gjg2fc+j _ J_ 

2n^^^* 2k 

j=0 1=0 j=o 


mrie+J{C) - log(n'^-C) ^ 


Define E 2 k-i{C) ■= J2^o which we note again is analytic at zero. Finally, 

defining c := — we see that by requiring consistency with 

(j4.32l) we have c = ne^. Let us also remark that using the definition of e together 
with {nih + '^^2 4“ 1 )^ 3 :* = 1 we see that ne^ = 0{1) as n —>• oo. 

To make statements concerning the properties of Ei near zero recall that by (11.2611 
we have. 


y+{x) G iM"*", if X G 5, (4.53) 

y{x) G M, if X ^ 5. 
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Note that the above implies hi{z), Ri{z) and E{z) are real functions for z G 
Consider, 

y(x* + C) = C)E{x^ + () (4.54) 

for X* + C in a compact subset M of \ {x*}. Note that E{x^) = 1 implies, for 
sufficiently small Dx,-, that E[x^ + C) ^ for C ^ MnM, hence we Hnd for sufficiently 
large n that, 

yi+(n‘^^*C) G *+, if X* + C G 5, (4.55) 

yi(n^^*C)GK, ifx*+C^5. 

Letting ( G M (1 M“ and using the large C expressions for yi in each case, gives, 

• (i) when x* is an interior point, 

CX) 

C)"' e (4.56) 

1=0 

which implies go € ^1^"*" and qi G iM. 

• (ii) when x* is a right edge point, 

CX) 

g«(A:+i)|^|fc+i ^g;(n^-*C)"' G (4.57) 

1=0 

which implies for k G that go G M"*" and g; G M, whereas for A: = — 1 it implies 
go G M“ and qi G M. 

• (hi) when x* is an exterior point, 

CX) 

1=0 

which implies g; G M. 

Finally, using the definition of Ei in terms of % and combining everything gives (j4.29p . 

□ 


4.3 Construction of 

Lemma 4.7. The local parametrix can be written in the form, 

P("*)(2) = P(z-x*)n'^"3$(/(z)|/,P,{rfe},foo)e™3®(^fo-, (4.59) 

where ^ is a canonical model problem with open lenses, E is an analytic function at zero, 
f is defined in Definition \1.21\ and G M. The value of the n-dependent parameters I, 
B, {rfo} and Too are given in the proof below. 

Proof. For clarity we split the construction into a number of distinct steps. 
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Matching the jump contours and singularities. 

First note that if <1> has a jump on a contour C then will have a jump on f~^{C). 

We therefore choose I = f{T\xt) and the lens contours to be such that = /(S±) 
for z £ Dx,. Similarly will have singularities at f~^{b) for b ^ B. To match 

the behaviour of T we therefore must choose B = f{B\x,) and for 

j > 0. This choice means that have behaviour 0(1) as n —)• oo due to Assumption 
a Note also that this choice is consistent with P^^*\z) having jumps on T;, for 
b £ B\x, since it can be verihed that /“^(Tb) = T/-i(6)- 


Matching the jump matrices 

We now consider the jumps of P^^*\z). By direct computation we hnd the jumps of 
(|4.59l) will match those of T if, for z £ Dx ^, we have, 


«r(/(^)) 


ar{z)-a+^ zGRHDx,, 

ariz) z £ nUfeeBU. 


By Definition 13.41 we have that for z £ Dx^, 


(4.60) 


^ j + EbeBU* Xr,{z)ab, zeRDDx, 

[EfeeBb* Xrtiz)ab, z £ Dx^ n UbeBU* \ 1^, 

and hence, 


(4.61) 


^ ^ fEbeBU.Xr,(/ ^elR, 

lEbeBU* Xr,{f~^{z))ab z £ IJbeBU. \ M. 


(4.62) 


Using the fact that Xri,{f~^{z)) = Xrf(^t,){z) and changing the summation dummy vari¬ 
able to 6 £ f{B\x^) we find that the jumps match if we choose dfe = aj-i^b) for b G B. 
We now have that (j4.59p fulfils conditions (a), (b) and (c) of the RH problem for 


Matching the boundary conditions on dDx, 

We now turn our attention to condition (d) of the RH problem for pP*'). We accomplish 
this by use of lemmas 14.31 and 14.51 We now consider each case in turn: 

(i) If x* £ 5*, i.e. an interior point, then for z = x* -|- C £ dDx^ we have, 

p(oo)(z)e’^?(^)'^3 ^ ('4^03) 

Taking to be of the form (I4.59p we have, 

= E{z)n^'^^x (4.64) 

= E(C)n^^^ {I + 0(n-^-)) QiC) (/ + ©(n"^-)) Q{C)-^E{C)-\ (4.65) 

= (/ + 0(n-^"*)) . (4.66) 

In the hrst line of the above we have used the asymptotics for <l>(z) as z ^ oo 
with Tooj = Ej{Q and roo,o = -®o(C) + (2fe + l)(m,^_|_(x*) -|- Tra^^). In the second 
line we have set P(C) = -®(C) and 4> = 0. 
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(ii) If X* G 55* with x* forming the right edge of an interval in 55*, i.e. a right edge 
point, then for 2 : = x* + C G we have, 

2fc+3 ESo(«^"‘C)"+^ *EdCV3_ (4_g7) 

Taking p(^*) to be of the form (|4.59l) we have, 

^ E{C,)n^^^ X (4.68) 

= E(C) (/ + C>(n-^)) (-f + iV-^C^^(C)”\ (4.69) 

= [l + 0{n-^)Y (4.70) 

In the first line of the above we have set (/> = A/4 and used the asymptotics for 
^>( 2 ) as 2 ; —)■ 00 with Tooj = EjiO- In the second line we have set E{Q = E{Q. 

(iii) If X* G 5* then, for 2 = x* + C G dDx^, we have, 

p(oo)(2)e"'^(^)°'3 = ii'(^)^/^'^3g7r*«7*^'(C)o'3+co-3log(n^^*C)-^ 

(4.71) 

Taking pP*) to be of the form ()4.59l) we have, 

(x* + ()e-<P*+f)'^3p{oo) ^ E{C)n^^^ X (4.72) 

(/ + n^“^**^“*°*+^“'^)‘^3^(atot+c-c-^)fT3g^Ei=2fe+i(»^^“^*C)^''“'.E;i(c)cr3 

= E{C)n^^^ (/ + 0(n-^"*)) (/+ ©(n"^"*)) P(C)"^ 

= ^/+ . (4.73) 

In the first line of the above we have used the asymptotics for d>(z) as 2 ^ 00 
with Tqoj{z) = Ej{z). In the second line we have used that /3 := — uef,^, 

dtot = «x*5 c = nef.^, — J\f and E{z) = E{z). In the third line we 

have chosen (f) = —'^a;*(dtot + c — c) = — A 3 ;^(Q!tot — 'nef;*). Note that in order for 
condition (d) to be satisfied in this case we require jatot — < 1/2. 

Let us note that the choices made for in each case are consistent with the existence 
of a solution to the model problem due to the properties, shown in Lemma 14.51 of the 
El in a small neighbourhood around zero. In particular we see that if is taken 
sufficiently small then Too will be in the neighbourhood of an admissible vector and 
therefore the RH is solvable (see Remark 12.81) and (I4.59P is well defined for our choice 
of Too. □ 

Remark 4.8. Note that we have again omitted the left edge case. The local parametrix 
at a left edge can be constructed in terms of the edge type model problem as follows, 

P(^*)(2) = P(2)ni^"*‘"3^3$(-/(2))a3e™30(P«7* _ (474) 

Checking the above satisfies the conditions of the local parametrix follows the case of 
X* being a right edge point. 
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5 Asymptotics for Kn{x,y): Proof of theorem 11.26 

We begin by expressing the correlation kernel in terms of the Y RH problem, 


Kn{x,y) 


1 ^yw{x)w{y) . 
27ri X — y ^ 


l)y+(y)-'F+(x) 



(5.1) 


The proof proceeds by inverting the sequence of transformations Y i—)-Ti—>-ri—)-i? 
in the steepest descent analysis. The result is that for 2 £ Dx^ we have, 

Y{z) = (5.2) 

Note that in the above equation we have used the fact that K= 

K{f{z)). By substituting the above equation into (|5.ip we arrive at. 


Kn{x,y) = -- 


gSTTioJ* 


I W\,r{x)whr{y) 


X 


(5.3) 


27ri(x -y)\j Wbr{x)whriy) 

Mf{y)))n-^‘^^E{y)-^R{yr^R{x)E{x)n‘^^^ {Mf{x)), Mf{x))f 

where we have introduced the functions 4)i{z) from Definition (ll.55p . Noting that. 


Wbrix) 
Wbr{x) 
we obtain, 


_ —27riar{x) 

cl • 


g7ri(2a+^ -ar(3;)-ar(y)) 

Kn{x,y) = -^- X 


(5.4) 


(5.5) 


27ri(x — y) 

{-<p2ifiy)), Mfiy)))n~'^‘'^E{y)-^R{y)-^R{x)E{x)n't’^^ {Mfix)), Mfix)))^ 


Next, using the properties of R, we have that, 

R(x* + n~^'^*v)~^R(x^ + n~^'^*u) = I + 0( ^ . ^ ). (5.6) 

Similarly we have that, 

E{x^ + n~^^*v)~^E{x^ + n~^^*u) = / + 0( ^. ^ ). (5.7) 

Using the above expression together with those for tcbr we obtain, 

Kn{x^+n-^-*U,X^+n-^-*v) = -arU*+n-^-*«)-arU*+n-^-*D) x 

0 l(tt) 02 (u) - (j)l{v)4>2{u) 

2'Ki{u — v) 

as n —)• oo. Finally we note, 

ar{x^ + n~^^*u) = ^ X6(Re (x* + n~^^*u))ab 
beB 

= ^ Xb(R-e{x^+n~^^*u))ab+ ^ Xb(R-e{x^ + n~^^*u))ab 

= ^ Xfe(Re (x* + n~^^*u))ab 

beBU* 

= + X] l^b(Re u)ab. (5.9) 

beB 


(1 + C>(n"^"*+^l‘^l)), 


(5.8) 
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We therefore obtain, 


n ^^*Kn{x^,+n 


— 


'u, X* + n = 


(5.10) 


^—iTi{a.Y ('u)+Q:r (i^)) 

27rz 


4>i{u)(I)2{v) - Mv)(I^2{u) + 0(n-^-+2|<^l)), 


u — V 


from which the theorem follows. The data for the model problem comes directly from 
the assignments made in the construction of the local parametrix. This is easy in the 
case of I, B and {Tb}. For Too note that from the construction of the local parametrix 
in the exterior and edge case we have that at finite n, Too,z = f{z)). In the 

case of an interior point we have Too,i{z) = f{z)) with the exception of roo,o 

for which we have, 


Too,o(^) = Eo{n ^"*f{z)) + {2k + l)(mC+(x*) + 7ra+). 


(5.11) 


The vector Too appearing in the model problem for the limiting kernel will be given by 
lim„^oowhich we now compute. 

For the edge and exterior case we have, 

lim Too(x* + n“^^*u) = lim Ei{0) (5.12) 

n—)-oo n^oo 

and for the interior case we have (j5.12l) for I ^ 0 and, 

lim Too,o{x* +n~^^*u) = lim (Flo(O) + {2k + l)(m,f+(x*) + vra^ )). (5.13) 

n—)-oo n^oo * 

Now note that, 

1(C) = lim yo(C) (5.14) 

n^oo 

which, using Ei{0) = qq^i, can be written for large C as follows, 

(i) if X* is an interior point of order k we have, 

2k+l 

= "2^ ^ ^,(0)C^' + 0(C-1), (5.15) 

i=o 


(ii) if X* is an edge point of order k we have, 

k-\-\ 

rk^(c) = + ^(C"'), (5.16) 

i=o 


(iii) if X* is an exterior point of order k we have. 


1(C) 


-clogC 


2k 


2k 

E 

1=0 


^,(0)C^+O(c-'). 


(5.17) 


This completes the final part of the proof showing that the model problem in Definition 
11.231 is the one appearing in Theorem 11.261 □ 
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6 Example applications 


6.1 The Painleve II kernel in the bulk 


We consider the model, 

( 6 . 1 ) 

■^n 

where V{M) = This model possess two distinct phases; for t > —2 the 

equilibrium measure support is connected, whereas for t < —2 the support consists of 
two disjoint intervals. The phase transition between these two phases was hrst analysed 
in [ 2 ^, in which it was shown that as the two cuts meet at the origin, the kernel there 
may be expressed in terms of Painleve II transcendents. We now will see how this result 
follows from Theorem 11.261 

From |26j we have the spectral curve, 

y{x) = {x^ + 2c) y/ — o? (6.2) 

with a = (|(—2t + 2{t^ + 12 ) 2))2 and c = |(t + + 8 ) 2 ). Consider the spectral 

curve when t = —2, 

y{x) = x^ — 4. (6.3) 


Now let t ——2 as n ^ 00 . In this limit we have 5* = [—2,2] and T-L^ = {0}. Using 
Defintion 1 1.1 71 we therefore see that 77* = {0}. So x* = 0 is in P* and from Definition 
11.131 is an interior point of order A: = 1. From Definition 11.181 we obtain For 

this model we have X = M and B = 0 which implies, using Definition 11.231 I = M and 
B = 0. 

_2 1 _2 

Now consider the limit n —)• 00 with t = —2 + n sr, this gives c = jTti 3 and we 
note that this means scales appropriately. Using Definition 11.221 we therefore have, 

y(C)=f(2C2 + r) (6.4) 


and 

l(C) = -f(^C' + rC). (6.5) 

Hence, from Dehnition ll.23l we have Too ,3 = 1, Too ,2 = 0 and Too,i = 3r. This completes 
the construction of the data for the canonical model problem. Note that this model 
problem is related, up to a rescaling of z, to the standard model problem for PII by 

(|2.18|) . Indeed we have, T^^^(cz) = e ^^(- 2 ) for ^ in the upper half plane and some 
_ *’^00,0 

constant c. Note that the factor e 3 cancels out in the expression for the kernel. 
From Theorem 11.261 we have that there exists a constant c such that. 


lim cn 3 Kn 

n^oo 


(n 3cu,n 



4'liu)(p2{v) - (Pl{v)4>2{u) 

2TTi{u — v) 


( 6 . 6 ) 


where. 


\Mz)) 


:= T^^^(z) 



(6.7) 


for z outside the lens. 
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1 0 
g27r7:(ai+a2) ^ 



Figure 8; The jump contours for the RH problem from which the new kernel at the 
hard-edge is constructed. 


6.2 A new kernel at a hard edge 

Consider the model, 

— I detMp“i| det(M - dM, (6.8) 

Zn 

where M is a negative definite hermitian matrix. We are going to consider the regime 
close to the origin as we allow t —)• 0 as n —)• oo. The matrix is taken to be negative 
definite so that the hard edge is a right edge point, which simplifies the construction. 

The spectral curve is the standard Marchenko-Pastur curve associated with the 
Laguerre unitary ensemble, 


y{x) = ( 6 - 9 ) 

From the expression for y we have 5* = [—4,0] and "H* = {0}. Using Defintion ll.171 we 
see that 7 ^* = {0} and therefore that the origin is a hard edge point with k = —1. From 
Definition 1 1.181 we obtain = 2. For this model we have X = and B = {0, t} with 
ao = ai and at = a 2 - We require t scales appropriately and therefore we set t = n~‘^T. 
This implies we have I = and B = {0, r}. The scaling limit of y is, 

y(c) = -2ru (6.10) 

and therefore, 

i = 2CX ( 6 . 11 ) 

which gives Tqo.o = — 1- 


This completes the construction of the data for the canonical model problem. To 
summarise, the model problem that will appear in the kernel is, 


44 









RH for 

(a) ; C \ M“ \ S* \ S® —)• is analytic in z. 

(b) $ has the jumps shown in Figure [8] 

(c) As z —)• oo, 

$(z) = (/ + Oiz-^)) 


(d) As 2 ; —>■ 0, 

$(z) = 0{l){z-h)^^^k{z). 


This can be written as, 


${z) = oi^^ 

\z 2 Z 2 


for z outside the lens and 


\z 2 z 2 


for z inside the lens. 
As Z ^ T, 


$(z) = O 


y(z-r)"^ 


{z-t)-^\ 

(z-t)-^J- 


( 6 . 12 ) 


(6.13) 


(6.14) 


(6.15) 


(6.16) 


Let us remark that this RH problem leads to Lax pairs which produce a second order 
ODE and therefore will be one of the standard Painleve equations. 

From Theorem 11.261 we have, 


lim n 3 Kn 

n—^oo 


3u,n 



where. 


[(f>2(z)J 



(t>l(u)(t>2(v) - (^l(v)(/>2(u) 

27ii(u — v) 


(6.17) 


(6.18) 


for z eM\I\Ub(=Brb- 
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